Markovianity criteria for quantum evolution 

Dariusz Chrusciriski and Andrzej Kossakowski 
Institute of Physics, Nicolaus Copernicus University, 
Grudziadzka 5/7, 87-100 Toruri, Poland 



Abstract 

We characterize a class of Markovian dynamics using the concept of divisible dynamical map. 
Moreover we provide a family of criteria which can distinguish Markovian and non-Markovian 
dynamics. These Markovianity criteria are based on a simple observation that Markovian dy- 
namics implies monotonic behavior of several well known quantities like distinguishability of 
states, fidelity, relative entropy and genuine entanglement measures. 

1 Introduction 

The dynamics of open quantum systems attracts nowadays increasing attention [H [21 [3]. It is 
relevant not only for the better understanding of quantum theory but it is fundamental in various 
modern applications of quantum mechanics. Since the system-environment interaction causes dis- 
sipation, decay and decoherence it is clear that dynamic of open systems is fundamental in modern 
quantum technologies, such as quantum communication, cryptography and computation [4j. 

The usual approach to the dynamics of an open quantum system consists in applying an ap- 
propriate Markovian approximation, that leads to the following local master equation 

Pt = Lpt , (1) 

where pt is the density matrix of the system investigated and L the time-independent generator 
of the dynamical semigroup. It is well known that under certain conditions on L [51 [6] the cor- 
responding dynamics p — > pt '■= Atp gives rise to completely positive and trace preserving map At 
[31 [1] (one usually calls a dynamical map) . The characteristic feature of Markovian approxima- 
tion leading to dynamical semigroup At = e tL is that it neglects all memory effects caused by the 
^ ■ nontrivial interaction of the system with an external world. However, recent theoretical studies 
and technological progress call for more refine approach based on non-Markovian evolution. 

Non-Markovian systems appear in many branches of physics, such as quantum optics [H [7], 
solid state physics [SJ, quantum chemistry [9], and quantum information processing [TO]. Since 
non-Markovian dynamics modifies monotonic decay of quantum coherence it turns out that when 
applied to composite systems it may protect quantum entanglement for longer time than standard 
Markovian evolution [llj. In particular it may protect the system against the sudden death of 
entanglement [12]. It is therefore not surprising that non-Markovian dynamics was intensively 
studied during last years jT3]— |2Tj . 

The standard approach to the dynamics of open system uses the Nakajima-Zwanzig projection 
operator technique [22] which shows that under fairly general conditions, the master equation for 
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the reduced system density matrix pt takes the form of the following non-local equation 

Pt = K-t-r Pt dr , 
Jo 



(2) 



in which quantum memory effects are taken into account through the introduction of the memory 
kernel /Q: this simply means that the rate of change of the state pt at time t depends on its 
history (starting at t = 0). It should be stresses that there is an perfectly equivalent approach 
[usually called time convolutionless (TCL)] [H [23] which describes quantum dynamics by local in 
time equation 



where Lt denotes local generator. It is clear that if Lt = L does not depend on time and ICt = 25(t)L, 
then both ([2]) and © reduce to the standard master equation ([1]). These equations provide, 
therefore, natural generalization of the standard Markovian semigroup. One of the open problems 
in this theory is to characterize properties of /Q and Lt which guarantee that the corresponding 
solution p t = A t p defines legitimate dynamical map A t . 

Interestingly, the concept of (non)Markovianity is not uniquely defined. One approach is based 
on the idea of the composition law which is essentially equivalent to the idea of divisibility |29j. 
This approach was used recently [31 1 to construct the corresponding measure of non-Markovianity. 
A different approach is presented in [30] where non-Markovian dynamics corresponds to a time 
evolution for the open system characterized by a temporary flow of information from the environ- 
ment back into the system and manifests itself as an increase in the distinguishability of pairs of 
evolving quantum states. The aim of this paper is to characterize a class of Markovian dynamics 
using the concept of divisible dynamical map and to provide a family of criteria which can distin- 
guish Markovian and non-Markovian dynamics, i.e. these criteria are satisfied if the dynamics is 
Markovian and hence their violation is a clear sign of non-Markovianity. 

The paper is organized as follows: in the next Section we recall the standard reduced dynamics 
of an open system. We stress that the problem of characterizing the properties of ICt and/or L t is 
in general untractable. It considerably simplifies in the case of commutative dynamics, i.e if the 
dynamical map At commutes in different times [A(,A U ] = for arbitrary t,u > (see Section [3]). 
In Section [J] we characterize Markovian dynamics using the concept of divisible dynamical map. 
Interestingly, we provide a simple method which enables one to go beyond Markovian dynamics 
fully controlling the corresponding local generator. This method is illustrated by pure decoherence. 
Necessary criteria for Markovianity are discussed in Section [6J Final conclusions are collected in 
Section [7J 

2 Reduced dynamics of an open system 

Consider an A^-level quantum system living in % coupled to a reservoir with the corresponding 
(usually infinite dimensional) Hilbert space Throughout the paper we use standard notation: 
B(H) and T(H) denote the class of bounded and trace class operataors in 7i, respectively. Let H 
denotes the Hamiltonian of the total composed system living in % (g> and uj be a fixed state of 
the reservoir. One defines the reduced dynamics At : T{%) — >T(%) by the following formula 



Pt = L t p t 



(3) 




(4) 
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where Ttr denotes the partial trace over the reservoir degrees of freedom. Note, that At is com- 
pletely positive and trace preserving for all t > and it satisfies Ao = 1. Therefore, it provides 
a legitimate quantum evolution of the system living in %. Actually, it is well known that any 
legitimate At may be defined as a reduced dynamics for appropriate Tin and the total Hamilto- 
nian H. The standard Nakajima-Zwanzig projection operator technique [221 [1] shows that reduced 
dynamics At satisfies the following non-local equation 

At = / Kt-T A rdr , A = 11 . (5) 
Jo 

The memory kernel Kit encodes all dynamical properties of the system and depends upon the total 
Hamiltonian of the "system + reservoir" and the reservoir reference state This equation is 
exact but in general very difficult to analyze. This is due to the fact that the memory kernel Kt 
depends upon all reservoir correlation functions. To simplify analysis one usually tries to perform 
suitable Markovian approximation to neglect all unwanted memory effects. The validity of such 
approximation is based on the existence of two characteristic time scales: the characteristic time 
Ts of variation of pt and the decay time tr of the reservoir correlation functions. Markovian 
approximation assumes that t$ S> tr. Basically there are two ways of rigorous treatment of the 
limit ts/tr —t- oo. One assumes that is invariant under the free evolution of the reservoir. 
Representing the Hamiltonian H as 

H = tf s ®I R + Is^#R + Atf int , (6) 

one performs the weak coupling limit A — > with rescaled time r = \ 2 t. In this scheme tr remains 
constant, while ts — > oo. This approach was analyzed in great details by Davies [241 125] (see also 
[26]). 

On the other hand in the singular coupling limit one has tr — > 0. It is achieved by considering 
the following Hamiltonian 

H = H s I R + e-% ®H R + e~ l H mt , (7) 

and performing the limit e — ¥ [27] . As a result reservoir correlation functions become 5 function. 
In both scenarios the limiting dynamics is governed by the well known master equation 

At = LA t , A = 1 , (8) 

where 1 denotes the identity map, and the Markovian generator L is given by 

Lp=-i[H,p] + ^([V a ,pV^ + [V a p,V^]) . (9) 

a 

In what follows we call L represented by ([9]) GKSL generator. In the above formula H represents 
the effective system Hamiltonian and {V^} is the collection of arbitrary operators encoding the 
interaction between system and the environment. Equation ([8]) gives rise to Markovian semigroup 
At = e Lt satisfying the following homogeneous composition law 

AfA u = At +n , (10) 
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for any t, u > 0. 

It should be stressed that one obtains Markovian master equation ([8]) from the general Nakajima- 
Zwanzig equation ([5]) only if the total Hamiltonian H enables one to perform suitable Markovian 
approximation and hence it covers only limited number of physically interesting systems. In general 
Markovian approximation is not suitable and one has to deal with much more involved non-local 
equation ([5]) . One of the main problems is to characterize the properties of the corresponding mem- 
ory kernel JCt which guaranties that the corresponding solution At represents legitimate dynamical 
map. 

Note, that instead of non-local equation ([5]) one may equivalently describes dynamics using 
local equation. This approach (usually called time-convolusionless pQ) leads to the following local 
in time equation 

At = L t A t , A = 1 , (11) 

with time-dependent local generator Lf We stress that these two approaches are equivalent. As- 
suming that At is differentiable it always satisfies local in time master equation. Indeed, formally 
one has At = A t A^ 1 A t = L t A t , where we assumed the existence of the inverse map A f _1 . Note, 
that the inverse, even if it exists, needs not be completely positive. Again, one would like to per- 
form characterization of time-dependent generators Lf giving rise to legitimate dynamical maps At. 
This problem seems to be untractable in full generality. Note, that formal solution to (jlip has the 
following form 

At = T exp (J L u di?J , (12) 

where 'T' stands for chronological product. It is clear that the above expression has rather a formal 
meaning. 



3 Commutative class of dynamical maps 

As we already stressed the general solution to the local in time master equation (lllj) has only a 
formal meaning and in general we do not control the properties of Lf which guarantee that T- 
product exponential formula Texp(J Q i L u du) defines dynamical map. Note, however, that if Lt 
defines a commutative family of generators, that is 

[L t ,L u ] = 0, (13) 

for any t, u > 0, the formula (112p considerably simplifies: the 'T' product drops out and the solution 
is fully controlled by the integral J Q * L u du. 

Theorem 3.1 If Lt defines a commutative family, then Lt is a legitimate generator of a quantum 
dynamical map if and only if J Q * L u du defines a legitimate GKSL generator for all t > 0. 

Example 3.1 As an example of commutative dynamics consider the following evolution of the 
qubit 

L tP = —^[<Jz,p] + ^-(cr z pa z - p) , (14) 
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where u>t, Jt '■ ^ + — ^K. The corresponding density matrix evolves as follows 

n-( P0 ° Poie [ ~ iQt ~ rt) \ n ^ 

Pt ~ I Pioe^ t) pu ) , (15) 

where ^ 

n t = / , r t = / 7„(iu . 

JO Jo 

L t gives rise to legitimate quantum evolution if and only if T t > for all t > and hence the 
evolution corresponds to simple decoherence. The corresponding dynamics is Markovian iff jt > 0. 

This example may be easily generalized for d-level system. Consider the following class of 
generators: let A = e 2m ^ d and define 

d-l 

V a = ^2\ a ?Pp, (16) 

p=o 

where Pa = |e^)(e^| and {eo, • ■ ■ , e^-i} stands for an arbitrary orthonormal basis in C d . Now, let 
Capit) be a time-dependent hermitian matrix and define 

d-l 

L t p = -i[H t ,p] + Cki(t) ([V k ,pV?] + \V k p,vfy , (17) 

k,l=l 

where the time-dependent Hamiltonian Ht reads as follows 

d-l 



Ht = + h k {t)V ] k ) , 

k=l 

where h k : M + —>C Lt generates legitimate dynamics if and only if the matrix 

Cki(t) = / c H (u)du 
Jo 

is positive definite. Dynamics is Markovian if and only if c k i{t) is itself positive definite. 

4 Markovian dynamics 

In this section we characterize important class of dynamical maps representing Markovian evolution. 
It is important to clarify this issues since there are few definitions used in the literature recently. 
We call a dynamical map At divisible if for any t > s > one has the following decomposition 

A t = V t>s A s , (18) 

with completely positive propagator Vt >s . Note, that Vt tS satisfies inhomogeneous composition law 

V t>s V SiU = V t , u , (19) 

for any t > s > u. In this paper following [31] we accept the following 
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Definition 4.1 Dynamical map At corresponds to Markovian evolution if and only if it is divisible. 

Interestingly, the property of being Markovian (or divisible) is fully characterized in terms of the 
local generator L t - Note, that if A t satisfies (fTT]) then V t:S satisfies 

d t V t , s = L t V t , s , V S;S = 1 , (20) 

and the corresponding solution reads 



ft 

V t a = T exp 



j Ludu^J . (21) 



The central result consists in the following 

Theorem 4.1 The map At is divisible if and only if Lt has the GKSL form for all t. 

Proof: (the proof goes similarly as for time-independent case (cf. [HE]) Assume that At is divisible, 
that is, Vt tS satisfies (|20p . One has 

Lt = lim V *+* ~ 1 , (22) 

for any t. Now, to compute Lt let F a (a = 0, 1, . . . , d 2 — 1) denotes an orthonormal basis in Md(C) 
such that F = I^/y/d. Now, since Vj jS is completely positive one has the corresponding Kraus 
representation 

d 2 -i 

V t>sP = c aP (t,8)F aP F} , (23) 

a,/3=0 

where the matrix c a p{t, s) is positive definite for all t > s , and V^t = 1 implies c a p{t } t) = dS a oSso 
for all a, /3 = 0, 1, . . . , d 2 — 1 . One finds 

( d 2 -l 

I coo(t + e,t) - d Tr-y cp k (t + e,t) t 
Lt p = lim < F pF + > F oP F k 



k=l 



E ?2 ^ ±M + E ^ > • (24) 



k=l k,l=l 



Let us introduce 



u\ y c 00 (t + e,t) - d c ok (t + e,t) c H (t + e,t) 

a 00 {t) = lim , a ok (t) = lim , a M (t) = hm ■ 



e^o e <=->o e e^o e 

Note, that the Hermitian time-dependent matrix au{t) is positive definite. Moreover, let 

Gt = Yd a ° o{t)Id + \ {At + A ^ ' Ht = Yi {At ~ A ^ ' 

with A t defined by 

At = —= E a k(t)F k . 
^ d k=i 
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Finally, one obtains the following formula for the local generator 



d 2 -i 

L tP = -i[H t , p] + {G t ,p} + a ki(t) F kP F} . (25) 

k,l=l 



Taking into account that At is trace preserving one has Ti(Ltp) = for all p which implies 

2 



Gt = -^a H (t)F^ , (26) 



k,l=l 
and hence 

d 2 -i 



L t p = -i[H t , p]+Yl [ F kP F i ~ \i F l F k> P}j . (27) 



k,l=l 



reproduces the standard GKSL form of Lt (recall that a,ki{t) is positive definite). 
Assume now, that Lt is defined by (|27p . It may be rewritten as follows 



Lt = $ t -* t , (28) 

where &t is a family of completely positive maps 

d 2 -i 

<5>tP = a ki(t)F kP F} , (29) 
k,l=l 

and 

%p = C tP - P C\ , (30) 

where 

C t = iH t + G t . (31) 

Actually, due to ([26]) one has Gt = — Note, that by construction Tr(L t p) = 0, and hence the 
corresponding solution Vt jS is trace preserving. It remains to show that Vf s is completely positive 
for all t > s. Consider the following equation 

d t N t , s = -^tNt,s , N s>s = t. (32) 

One easily finds 

N tjS p = X ttS pXl, (33) 
where X t)S itself satisfies dtXt jS = CtX tjS and hence 

Xt, t = Texp(J C u di?J . (34) 

It is clear that Nt 8 is completely positive. Moreover, it is invertible and the inverse Nf s reads 

Nr}p = Y ti ,pYl , (35) 
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where 

r-t 



Y t>a = T exp I - J C u duj , (36) 

with To denoting anti-chronological operator. Hence is completely positive as well. To solve 
original equation ([20]) let us pass to the "interaction" picture and define 

V t>3 = NttVg? . (37) 

One finds 

dtV ^ = ^V^, V%0 = 1, (38) 



where 

- t,s 



C^^^i^M. (39) 



The above formula shows that $?f " is completely positive being the composition of three completely 
positive maps: N^ 8 , $>t and Nf s . One easily solves ([38]) and gets 

= T exp (Y du\ . (40) 
It is therefore clear that can be represented as the following series 

^i nt) = t+f d h + f dtl T dt2 o *gy + . . . . (4i) 

It shows that V t s is completely positive being a sum of completely positive maps. Hence, taking 



into account formula (j37j) it finally shows that Vt tS is completely positive. □ 
It 

v} s . Note, that formula (|39j) implies 



It is clear from (|4ip that complete positivity of is sufficient for complete positivity of 



^ = N*o^f*oN*-\ (42) 

where A# : B{T-L) — >■ B{T-L) denotes a dual map defined by 

Tr(p • A # a) = Tr(a • Ap) , (43) 

for all p £ T(J-L) and a G £>(%). Now, if <I>[ mt ^ is completely positive then formula (I42p proves that 
<J>i is completely positive as well. Hence, complete positivity of <E>i cannot be relaxed. 

Let us observe, that a class of Markovian evolution may be easily generalized as follows: we 
are looking for the solution of (jlip with Lt represented as in (|28p . Our aim is find <I?t and ^ 
such that L t defines a legitimate generator. Suppose, that we are given a completely positive map 
Nt satisfying No = 1. Suppose that Nt is not trace preserving (if it were it represents legitimate 
dynamics and we are done). Let us define *&t as a local in time generator such that 

N t = -%N t , N = t. (44) 
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One obviously has ^ = —N t N t 1 . Now, the question is: does there exit $ t such that L t = <& t — ^ t 

(int) 

4 > 



provides legitimate generator? Note, that the solution of (jlip has the following form A t = N t A^ mt 



where A| mt ^ is defined in terms of the following series 



A? nt) = i + f dh d>J; nt) + f dtx I " l dt 2 < nt) o 4^ + . . . , (45) 

Jo Jo Jo 

with $| int ^ = Nf 1 o $ t o Nt- It is clear that if <3>£ int ^ is completely positive so is Aj. Now, 

= JVt o $f nt) o A^ 1 = QtNf 1 , (46) 

with @t := At$j being a completely positive map. Note, that <&t needs not be completely 
positive. However, if Nf is completely positive then necessarily <3?t is completely positive as well. 



Hence, if $t is constructed by (j46J), then Af is completely positive. It remain to check for trace 
preservation. Note, that At is trace preserving if L^I = 0. One has 

L*l = N t - 1 *(e* + N*)l, (47) 

and hence if 

Of I + N*I = , (48) 
then At defines legitimate dynamical map. It is therefore clear that if 

- N*I > , (49) 



then one can always find completely positive ©t such that normalization condition (I48D holds. 
Clearly, the choice of ©t is highly non unique. If ©f satisfies (|48p and Mt is an arbitrary family 
of quantum channels, then the following 'gauge transformation' @t — yQf 1 := Mt@t gives rise to 
another admissible Qf 1 satisfying ([4"Hj). 

Proposition 4.1 If Nt with Nq = i is a family of completely positive map satisfying ( f^ff[ ), then 
there exists completely positive ©t satisfying f^<g| ) such that Lt = (©t + N^Nf 1 gives rise to the 
legitimate local generator. 



Actually, condition (J49J) is always satisfied for the Markovian evolution. Indeed, taking into account 
([33]) one finds 

d t N*i = (d t xl)x t , s + xl(d t x t , s ) , (50) 

and hence using dtX t , s = CtX ttS one finds 



d t N#I = Xl 



Ct + Cf A tjS — 2Xj. s GtX tjS , (51) 



-t 



which shows that dtN^ s I < due to Gt < 0. Therefore, presented method generalizes Markovian 
generator keeping N t completely positive and satisfying (|49j) but admitting & t to be not com- 
pletely positive. It should be stressed that this construction provides a local analog of semi-Markov 
dynamics constructed recently in |34| . 
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Example 4.1 Consider the following family of completely positive maps 

N 

N t p= ^2 n ki(.t) ekk P en , (52) 
k,l=l 

where = and the matrix n k i(t) is positive definite with Wjy(O) = 1. One easily finds for *$>t 

N . , s 

ftp := -N t N^p = -J2 ^77T efcfc P en . (53) 
fc,z=i j 

Note, that condition \J^9ty is equivalent to the following condition for the diagonal elements 

h kk (t) < , (54) 
which implies n,kk(t) < nkk(0) = 1. Now, let @t be a family of completely positive maps 

N 

&tP= ^2 °ki{t)e kk peu , (55) 

k,l=l 

where the matrix 9ki{t) is positive definite. Normalization condition fl^gp shows that the diagonal 
elements of the matrix 0ki(t) are uniquely determined by Okk{t) = ~ n kk{t) ■ The off-diagonal ele- 
ments 6ki(t) are arbitrary provided that Okiif) is positive definite. The simplest choice corresponds 
to 9ki{t) = for k I. It finally gives 

N . ,s 

$ t p = e t Nf l p = V e kk P e kk , (56) 

l~y n kk{t) 

and hence 

L t p=(<S> t -%)p = y2 e kk P e u , (57) 

provides the pure decoherence dynamics: A t p = J2ki=i c ki{t) ekk peu with 

Cki{t) = n k i(t) , (k^l) and c k k = 1 • (58) 
The matrix Cki(t) is by construction positive definite. 



5 Characterizing Markovian dynamics 

In this section we analyze special properties of divisible (and hence Markovian) dynamical maps. 
Let us recall, that if a linear map A : T(T~L) — > T(7i) is trace preserving, then A is positive if and 
only if 

||Aa||i < ||a||i , (59) 

for all hermitian a. Note, that A needs not be contractive for non-hermitian elements. However, if 
A is completely positive, then (|59p holds for all a G B(H). Actually, it turns out [35] that if A is 
2-positive and trace preserving, then 

||A||i := sup ||Ao||i = sup ||Aa||i = 1 . (60) 

||a|ji=l ||a||i=l;a+=a 
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Corollary 5.1 If At is a dynamical map, then ||Aj||i = 1, that is, At is contractive in the trace 
norm 

IIAHIl < IHIi ■ (61) 

Moreover, if At is a divisible map, then 

^l|A 4 a||i<0, (62) 

for an arbitrary a £ T(rl). 

Similar property holds for dynamical maps in the Heisenberg picture. Recall, that if A t is a 
dynamical map in the Schrodinger picture, then its dual A^ : B(H) — > B(H) corresponds to the 
Heisenberg picture. It is clear that A^ is a unital completely positive map for all t > 0. If A* is 
unital and completely positive, then its operator norm satisfies 

||A # || := sup ||A # a|| = 1 . (63) 
IMI=i 

Corollary 5.2 If Af is a dynamical map in the Heisenberg picture, then ||A*|| = 1, that is, A* 
is contractive in the operator norm 

\\A*a\ \ <\\a\\ . (64) 



Moreover, if A t is a divisible map, then 



|||Afa||<0, (65) 



for an arbitrary a G B(rL). 

Example 5.1 Consider once more the generator defined in 1 1J$ . One has 



where cr+ = |0)(1|. It implies 



A*a+\\ = \\e int - Tt a + \\ = e- rt ||cj + || = e~ Tt , (66) 



^||Afa + ||=-f t = - 7 tr* , (67) 



which shows that Markovianity of At implies "ft > 0. 

Let us observe that if the total Hamiltonian H of the "system + reservoir' has a discrete spectrum 

H = ^e a P a , (68) 



then dynamical map At defined by 

A t p = Tr R 

has the following form 



(p®tJ R )e im , (69) 



A tP = y £e-^-^ t A a pp, (70) 

a,/3 
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where A a pp = Tr n(P a [p® u)R]Pp). It is clear that due to the presence of the oscillatory terms 
e -t{t a -ep)t corresponding trace norm ||A t a||i is an almost quasi-periodic function and hence can 
not be monotonically decreasing. It proves that in such a case one obtains genuine non-Markovian 
dynamics. 

If A : T(H) -^T(H) is a linear map one defines so called diamond norm 

||A|| := sup ||(l®A)W||i . (71) 
||W||i=l 

Theorem 5.1 Let A t be a dynamical map. The following conditions are equivalent 

1. At is divisible, 

2- |H4, a ||o = 1 for all t > s, 

3. one has 

^IKlO^Wlli^O, (72) 

for all Hermitian W € T{H®%). 

The corresponding theorem in the Heisenberg picture may be formulated as follows: recall that 
A* : B(H) -+B{U) is completely bounded if 

||A#|| cb :=||l®A#||<oc, (73) 

and it is completely contractive if ||A#|| c b < 1. 

Theorem 5.2 Let Af be a dynamical map in the Heisenberg picture. The following conditions are 
equivalent 

1. Af is divisible, 

2. ||V t *|| c b = 1 for all t> s, and hence V^ s is completely contractive, 

3. one has 

|||(l®Af)A||<0, (74) 

for allAe B(H®H). 

The Markovian evolution may be characterized in a slightly different way: we know that At corre- 
sponds to Markovian evolution iff the 2-parameter family of propagators V^ s is completely positive 
for t > s. Denote by ip + maximally entangled state m.7-L®l-L and let P + = Note that 

Vt tS is completely positive if and only if (1(g) Vt, s )P + > which is equivalent to the following simple 
condition 

v t , s := 11(1®^)^ ||i = 1 • (75) 
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Let us define 

dv ut 

9t ■= —j— 

du u=t 

that is, 



(76) 



9t = lim V J±^1 , (77) 
where we have used vtt = 1- Taking into account that 

F t+M = 1 + eLt + 0(e 2 ) , (78) 

one finds 

5* = I™ . 79 

e->0+ e 

Corollary 5.3 ( [31 J ) A map At is divisible if and only if gt = for all t > 0. 

Finally, let us provide characterization of the corresponding generator in the Heisenberg picture. 
Let us recall that if A* : B(H) — > B(H) is unital and 2-positive the following Kadison inequality 
holds 

A*{aa*)>A#{a)A#{a*). (80) 

This inequality may be used to characterize Markovian generators. Note, that Markovian dynamics 
A* satisfies 

d t V# = V#Lf, V* = t, (81) 
where V^ s denotes the dual propagator. Now, differentiating the Kadison inequality 

V*(aa*) > V*(a)V*(a*) , (82) 



one finds 

V#L*(aa*) > V*L#(a) ■ V*(a*) + V#(a) ■ V*L#(a*) . (83) 
Taking t = s and using V t t = 11 one gets 

L*(aa*) > L*(a) ■ a* + a ■ L*(a*) . (84) 
Definition 5.1 A hermitian map 9 : B(H) — > B(H) is dissipative iff 

9(aa*) > 9(a) ■ a* + a ■ 9(a*) , 
for all a G B(H). 9 is completely dissipative if ]l<g)\E f is dissipative. 
One has 

Proposition 5.1 A^ defines Markovian dynamics in the Heisenberg picture if and only if is 
completely dissipative. 
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6 Simple criteria for non-Markovianity 



In this section we develop a series of criteria for non-Mar kovian dynamics. It turns out that 
dynamics represented by a divisible map displays characteristic monotonic behavior for several 
interesting quantities. Breaking monotonicity reveals non-Markovian character of the corresponding 
quantum evolution. 



6.1 Distinguishability 

Trace norm defines a natural distance between quantum states represented by density operators: 
given two density operators p and a one defines 

DM = ^\\p-<r\\i. (85) 

The quantity D[p,a] is usually interpreted as a measure of distinguishability of quantum states p 
and a. It is well known that if A is a positive trace-preserving map, then 

D[Ap,Aa] < D[p,a] . (86) 
Corollary 6.1 If At is a divisible map, then 

j t D[A t p,A t a] <0 , (87) 

that is, for the Markovian evolution distinguishability of any pair of initial states monotonically 
decreases. 

It is well known that if At corresponds to the unitary dynamics Atp = UtpJjj , with unitary Ut, then 
D[Atp, Ato~] = 0. Moreover, if At = e Lt represents dynamical semigroup, then D[Atp, At<r] < 0. The 
above property was used by Breuer et al [30J as another definition of Markovianity. This criterion 
identifies non-Markovian dynamics with certain physical features of the system-reservoir interaction. 
They define non-Markovian dynamics as a time evolution for the open system characterized by 
a temporary flow of information from the environment back into the system. This backflow of 
information may manifest itself as an increase in the distinguishability of pairs of evolving quantum 
states. It turns out that these two concepts of Markovianity do not agree (see e.g. [331 E2])- Clearly, 
divisibility implies (i87|) but the converse needs not be true. 



6.2 Fidelity 

Given two density operators p and a one defines Uhlmann fidelity 



Equivalently, one has 

F(p,a) = \\^\\l, (89) 
which shows that F(p,a) = F(a,p). One proves 



1 - F(p, a) < D[p, a]<y/l- F(p, a) 2 . (90) 
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Moreover, If A is a quantum channel, then 

F(p,a)<F(Ap,Aa) 
Corollary 6.2 If At is a divisible map, then 



d 
dt 



F(A tP ,A t a) > 



6.3 Entropic quantities 

Let us recall the definition of Renyi S a and Tsallis T q relative entropies 



for a G [0, 1) U (1, oo), and 



Sol{p\\o) 



T q(p\\<r) 



1 



a — 1 
1 



log 



1 



l-Trp q a 1 - q 



for q G [0, 1). Note, that in the limit 



(91) 



(92) 



(93) 
(94) 



lim S a (p || a) = lim TJp\\a) = S(p \ \ a) , 

a — y 1 q — ¥ 1 

one recovers well known formula for relative entropy 

S(p\\a) = Tr(p [log p- logcr]) . 
It turns out [U [25] that if A is a quantum channel then S a and T q satisfy 

S a (Ap\\Aa) < S a (p\\a) , T q (Ap\\Aa) < T q (p\\a) , 
for a G [0, 1) U (1, 2] and q G [0, 1). Clearly, the same property holds for the relative entropy. 
Corollary 6.3 If At is a divisible map, then 



(95) 



(96) 



S a (Atp\\A t a) < 



j t T q {A tP \\A t a) < 



for a G [0, 1) U (1, 2], q G [0, 1), and 



6.4 Entanglement measures 



-S(A t p\\A t a)<0 . 
dt 



(97) 



(98) 



Consider a composed system living m.T-L®%' and let be an arbitrary density matrix in T-L^Ti'. 
It is well known [3] that for arbitrary genuine entanglement measure £ one has 



£([$®&}W) < £(W) , 



(99) 
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where $ : T(H) — >• T(H) and <J>' : T(H') — > T^H 1 ) are quantum channels. Denote by Wt the 
trajectory W% = (At (g) 1)W starting at W. Now, if E is an entanglement measure, then 

£{W t ) < £{W) , (100) 

for an arbitrary dynamical map At. It is, therefore, clear that if At is divisible then 

j t £(w t )<o, (101) 

for each initial state W. Of course the above relation is nontrivial only if the initial state W is 
entangled. In particular if W = P + (maximally entangled state), then entanglement measured by 
£ monotonically decreases from the maximal value £(P + ) |31j . 

7 Conclusions 

We characterized a class of Markovian dynamics using the concept of divisible dynamical map. 
Interestingly Markovian dynamics is fully controlled in the local approach via the properties of the 
corresponding local generator. Characterization of Markovianity in terms of the memory kernel tCt is 
an open problem. It should be stressed that the standard Markovian master equation ([I]) is defined 
by the corresponding macroscopic model and a suitable Markovian approximation. Note, that for 
general Markovian evolution characterized by time-dependent local generator the construction of 
the corresponding macroscopic model is not known. 

Moreover we provided a family of criteria which can distinguish Markovian and non-Markovian 
dynamics. These Markovianity criteria are based on a simple observation that Markovian dynamics 
implies monotonic behavior of several well known quantities like distinguishability of states, fidelity, 
relative entropy and genuine entanglement measures. 

We stress that the problem of characterization of admissible Lf and JCt giving rise to legitimate 
dynamical map At is rather untractable. Only commutative case is fully controlled. One may 
wonder if there is another way to describe dynamics of an open system. In a forthcoming paper we 
propose a new approach to this problem. 
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